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Recent theoretical studies of chaotic scattering have encounted ensembles of random 
matrices in which the eigenvalue probability density function contains a one-body 
factor with an exponent proportional to the number of eigenvalues. Two such 
ensembles have been encounted: an ensemble of unitary matrices specified by the 
so-called Poisson kernel, and the Laguerre ensemble of positive definite matrices. 
Here we consider various properties of these ensembles. Jack polynomial theory is 
used to prove a reproducing property of the Poisson kernel, and a certain unimodular 
mapping is used to demonstrate that the variance of a linear statistic is the same as 
in the Dyson circular ensemble. For the Laguerre ensemble, the scaled global density 
is calculated exactly for all even values of the parameter (3, while for (3 = 2 (random 
matrices with unitary symmetry), the neighbourhood of the smallest eigenvalue is 
shown to be in the soft edge universality class. 



1 Introduction 

In the theory of random matrices, a primary task is to compute the probability density function 
(p.d.f.) for the eigenvalues from knowledge of the p.d.f. for the ensemble of matrices. Two 
examples of random matrix ensembles of interest in this paper are Dyson's | Dys62| 1 circular 



ensembles of symmetric, unitary and self dual quaternion unitary random matrices (labelled by 
= 1, 2 and 4 respectively), and the Laguerre ensemble of random Wishart matrices A = X'X, 
where X is a random M x N (M > N) matrix which has either real {(3 = 1), complex ((3 = 2) 
or quaternion real {(3 = 4) Gaussian random elements. In the circular ensemble, the p.d.f. for 
the matrices is uniquely specified by requiring that it be uniform and unchanged by mappings 
of the form U i— > VUV', where V is an arbitrary N x N unitary matrix and V = V T for 
(3 = 1, V is arbitrary for (3 = 2 and V' = V D for (3 = 4 (D denotes the quaternion dual). 
In the Laguerre ensemble, the distribution of the elements of X are taken to be the Gaussian 
^ e -/3Tr(xtx)/2 -^high i s equivalent to choosing each element independently with a Gaussian 
distribution A'e~^ x ^ 2 / 2 . 

The corresponding p.d.f. for the eigenvalues e tSj , j = 1, . . . , N in the circular ensemble is 

J_ TJ | e ^_ e «Y ; (1.1) 

C(3N l<j<k<N 
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while for the Laguerre ensemble the p.d.f. is given by 

i N 

J-He-^xf 2 n \>*->*\ p , ^>°> a- 2 ) 

a0N j=i l<j<k<N 

where a := M — N + 1 — 2/(3. In both cases the eigenvalue p.d.f.'s can be written in the form 
of a Boltzmann factor for a classical gas, with potential energy consisting of one and two body 
terms only and interacting at inverse temperature /3, 

exp( -p^Vi( Xj )+ £ V 2 (x j ,x k ))\. (1.3) 

^ j=l l<j<k<N ' 

Thus for the circular ensemble 

V 1 (x)=0, V 2 (x,y) = -log\e ix -e% (1.4) 
while for the Laguerre ensemble 

X Ob 

Vl( x ) = 2 ~ 2 logX ' ^(a:,?/) = -log|x-y|, (1.5) 

which displays the well known fact that the analogous classical gas has a two body logarithmic 
potential. The logarithmic potential is special in that it is the Coulomb potential between like 
charges in two dimensions. In (|0|) the two dimensional charges are confined to a unit circle, 
while in (|0|) the two dimensional charges are free to move on the half line x > 0, but are 
confined to the neighbourhood of the origin by the one body potential. 

In this work we will focus attention on a subclass of random matrix ensembles with eigen- 
value p.d.f.'s of the form ( |1.3[ ) in which the one body potential contains a term proportional to 
iVlog |1 — fi*e tx \ (unitary matrices) or —Nlog \x\ (Laguerre ensemble). Again the two body term 
is logarithmic. In the classical gas these one body potentials can be interpreted as being due 
to an external fixed charge at the point 1//U* in the complex plane with strength proportional 
to —N (unitary matrices), and as an external charge fixed at the origin of stength proportional 
to N (Laguerre ensemble) . We see from (|1.2| ) that an example of a random matrix of this type 
in the Laguerre case is a Wishart matrix with X a rectangular matrix in which the number of 
columns is some fixed fraction of the number of rows. In the theory of random unitary matrices, 
this type of eigenvalue p.d.f. results as a special case of the the ensemble of random matrices 
defined by the Poisson kernel 

i det(i - gjtffy-iyw 

C\det(l-S^S)\^ N - l )+ 2 ' { ' 

where the notation S denotes the average of S. Random unitary matrices with this p.d.f. occur 
in the study of scattering problems in nuclear physics [MPS85] and mesoscopic systems |Bro95]. 



In the case 5 = 0, this reduces to the p.d.f. specifying Dyson's circular ensemble (all members 
equally probable). In the case that S = [iIn, |m| < 1 (In denotes the N x N unit matrix), the 
corresponding eigenvalue p.d.f. is given by 

l \ \1-I2*e i6 \^' , , v 

J=l 1 ^ 1 l<J<k<N 



C P n |l-/i*e^' 11 1 1 ( ' 
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where a' = (N — 1 + 2/(3) and Cpisr is as in and so in the classical gas picture there is 

a fixed charge of opposite sign at the point 1///*, and the magnitude of this charge is indeed 
proportional to N. 



In Section 2 we will consider the eigenvalue p.d.f. (|1.7| ) for the Poisson kernel. First the 
special reproducing property of the p.d.f.'s ( |1.6[ ) and ( |l~7| ) for (3 = 1,2 and 4 will be revised, 
and extended to all /3 > in the case of (|Oj) • Our tool here is Jack polynomial theory [ Mac95| ] . 
Then we will consider the effect of the iV-dependent exponent in (1.7) on the one and two body 
correlation functions, as well as on the fluctuation formula for the variance of a linear statistic. 

In Section 3 a physical problem giving rise to the Laguerre ensemble in which a in ( |1.2| ) is 
equal to N will be revised. Then we will revise known theorems in the cases (3 = 1 and 2 for the 
global density and the distribution of the largest and smallest eigenvalues. Next known integral 
formulas for the density [BF97a] at general even (3, deduced from the theory of generalized 
Selberg integrals [Kan93| and their relationship to Jack polynomial theory, will be analyzed in 
the appropriate limit to deduce that the formula for the global density holds independent of 
(3. For the special coupling (3 = 2 the local distribution functions in the neigbourhood of the 
smallest eigenvalue are analyzed and shown to belong to the universality class of the soft edge, 
giving rise to the Airy kernel [ |For95 , TW94a, KF97b]. We conclude the section with an analysis 
of some nonlinear equations [TW94a, TW94b], which explicitly demonstrates the universality 
of the distribution of the smallest eigenvalue. 



2 The Poisson kernel 

2.1 Physical origin of the Poisson kernel 

The scattering matrix for N channels entering and leaving a chaotic cavity via a non-ideal lead 
containing a tunnel barrier has as its p.d.f. the Poisson kernel (|1.6| ) [Bro95]. Also, in scattering 



problems in nuclear physics, ( |1.6| ) has been used |MPS85] to describe situations in which the 
average of the scattering matrix is non-zero. It was in the latter problem that the Poisson kernel 
first appeared in an application of random matrix theory. In | MPS85[| a requirement for the 



p.d.f. of the ensemble of scattering matrices was the special reproducing property 

f(S) := f((S)) = (f(S)) (2.1) 
for / analytic in S. It was noted that a result of Hua |Hua63[ gives that 

f(S) = cJ fW l<isb{ i N -s\s)\w-w KdS) (2 - 2) 

where fi(dS) is the invariant measure associated with the Dyson circular ensemble, and thus 
that the Poisson kernel exhibits the reproducing property fl2.ip. 
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2.2 The reproducing property for general j3 



To make any sense out of Q2.2| ) it is necessary that /? = 1,2 or 4, so that the measure /j,(dS) 
has meaning. However, in the case S = /iljv, the corresponding eigenvalue p.d.f. (|1.7|) can be 
interpreted as a Boltzmann factor and it makes sense to consider all > 0. In this case, for 
(3 = 1,2 and 4, ( |1,7| ) used in ( |2.2j ) gives that the reproducing property restricted to analytic 
functions of the eigenvalues e ,ej reads 



/(/v--,m) 



1 



N 



2- 



1 



|1 — /x*e 



„|2\/3a'/2 



* piQi I /3a' 



n 



(2.3) 



l<j<k<N 



This equation is well defined for all /3 > 0, but its validity has only been established for (5 = 1, 2 
and 4. Here we will establish its validity for general /3 > 0, using properties of the orthogonal 
polynomials associated with the p.d.f. Ql.ip . Note that in the case N = 1 (2^) reads 

H 2 ) 



i 

2^ 



277 



1 



(2.4) 



which is the celebrated Poisson formula on a circle, giving the value of an analytic function / for 
< 1 in terms of its value on the unit circle. Thus ( |2.3| ) can be regarded as an iV-dimensional 
generalization of this result. 

The multivariable orthogonal polynomials corresponding to ([O]) are the symmetric Jack 
polynomials P^ 2/l3 \z) jSta89| , |Mac95| , as they possess the orthogonality property 



N 

n 

i=i 



2tt 



l<j<k<N 



t8 k 



Jdj 1/3 _ Kf X 



(2.5) 



and form a complete set for the space of analytic functions. Here the labels k and a are partitions 



consisting of N parts or less, and z : = (e 

Mo 



£0i 



The normalization N K is given by 



[l + P(N-l)/2# 



(2//3) 



(2.6) 



where 



N 



< : =n(|«( s )+^)+i), Mi a) :=fi T[l! 



1) + K 4 



r(n-f(j-i) 



i=l - V" 2 

The notation s £ k refers to the diagram of the partition k, and a(s) = k. 



j is the arm length 



while l(s) 



i is the leg length (k' refers to the conjugate partition of k); see e.g. [Mac95]. 



The normalization Mq is the same quantity as the normalization Cp^ in ( |1.1| ), and has the 
explicit value 



f3N 



09/2)! 



(2.7) 



The quantity Pi 2// ^ (1^) in ( [2.6| ) also has an explicit evaluation [Sta89, Mac95 , BF97b], but it 
suits our purposes to leave it unevaluated. 



4 



With these preliminaries, we now pose the problem of specifying the kernel K(w,z*), w = 
(e*^ 1 ,. . . , e l< ^ N ), Im((j)j) > 0, which is an analytic function of z*, and has the reproducing property 

(2.8) 

l<3< 

for / analytic and symmetric in w\, . . . , iojv- The proof of ( |2,3| ) will then consist first of evaluating 
K(w, z*) at w = (fx, . . . , fj,), then transforming the equation ( pH| ) so that K(w, z*) becomes real. 
Using the orthogonality property fl2.5|) , as well as the completeness of {P^^\z)} for analytic 
functions, it is a simple exercise to check that the required kernel K is uniquely given by 



N ,277 

/H = II / dBiK(w,z*)f(z) J] \e l - 

l=l J ° l<j<k<N 



K(w,z* 



E 



pi 2/ ^HPi 2// V 



1 



53[l + / 3(iV_l)/2]P 



(2//3) 



P, 



(2.9) 



where the second equality follows from ( |2.6[ ). But in general the generalized hypergeometric 
function 1 (a; u;, z*) is defined by 

r(2//3) 



(2//3) f K [W)"k [Z 



so we have 



1 



K(w, z*) = -^ 2 //3) (1 + 0{N _ 1)/2; w . 



(2.10) 



(2.11) 



In general there is no known expression for^^a ; w, z*) in terms of elementary functions. 
However, the case w = (fi, . . . , fj,) is an exception, for then we have 

M N [O ]£/0 P MV) 



4 2m (a;w,z*) 



W=(fl,...,/J,) 



dL 



(2.12) 



The significance of this is that \Jq {p\ fiz*) can be summed according the generalized binomial 

1 



formula [Kan93] 

N 

Comparing ( gH|) - (PH) we therefore have 



#(10,2*) 



if=(M>->^) 



1 * 



1 



^1)1+^-1)72 • 



(2.13) 



(2.14) 



3=1 v ~ r J- 

Although this is an explicit solution to the problem of determining the kernel in (2.8) in the 
case w = (fi, . . . , ju), it does not immediately establish ( |2.3[ ) as the kernel (2.14) is not real. Note 
that in the case N = 1, ( 2.14 ) is the Cauchy kernel from elementary complex analysis. For 
general N, to obtain a real (Poisson) kernel from the Cauchy kernel, we proceed as in the one 
dimensional case and simply make the replacement 



N 



1 



r/ 



A ! (l-M%-) 1+/3(Ar_1)/2 ' 

in ( j2.8|) with w = (fi, . . . , /x). The formula ( |2.3|) results, thus establishing its validity for general 
>0. 
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2.3 Fluctuation formulas 



In the application of random matrix theory, an important class of observables are the linear 



statistics A = J2j=i a(Xj). The first two moments of these statistics are given by 



2k 



p {1) (9)a(9)d9 



2tt 



d6 ia { 



2- 



d9 2 a(9 2 )S(9 1 ,9 2 ) 



(2.15) 
(2.16) 



(A) = 
Vai(A) : = 

where S(Xi, X 2 ) denotes the structure function 

S{X 1 ,X 2 ) :=pf 2) (Ai,A 2 ) + P(i)(Ai)<J(Ai -A 2 ), 

with p(i) denoting the density and pj 2 ^ denoting the truncated two particle distribution function. 
In this subsection these quantities will be considered for the eigenvalue p.d.f. ([L 



It is instructive to first consider the case p = 0, when ( |l.7|) reduces to the eigenvalue 
p.d.f. (1.1) for the circular ensemble. In this case [For9E] 

P(i)(* (0) ) = N 



2tt 



1 



d 2 



^ 2 ae^deP 



log I sin( 



,(0) 



4°))/2| 



(2.17) 
(2.18) 



(the use of the superscript (0) indicates that p = in (1.7)), where in the asymptotic expression 
(|2.18| ) all oscillatory terms, each of which have period some integer multiple of 2tt/N, are ignored. 
For example, we have the exact result (see e.g. JMch91[ |) 

f 1 ^ sin 2 N(9 { ° } -9 { 2 0) )/2 

1 cos N(9 ( ° ) -9 { 2 



n T (() (o) fl (o) N 



V2vr/ 



2\2vr/ V 



sin 



3 (0) 



Ignoring the oscillatory term with the factor of cos N(9i >> — 9^) gives 



4°))/2 

j(0) 



sin 



P(2)\Vl > 



(0) fl (0)x 



2 



N^oc 2(2^)2 sin 2 ( ^(0) _ (o) )/2 " 2n 2 de f) de 



(0) 



log I sin( 



j(0) 



vaild for 9^ 7^ #2°' ) - The validity of ( |2.18 ) for 9^ = 62 is then deduced by noting that if we 
write (define) 

i-2-k Q2 



j(0) 



3 (0) 



log|sin(^ 0) -^ 0) )/2|^i 



oe¥>oof> 



3(0) 







d9 { 2 0) Jo d9\ 



2- 



d 1 I • fa(0) 
— log I sm(#i 



9 { 2 0) )/2\d9 ( "> (2.19) 



(0) 



(0) 



then this integral vanishes. This is the perfect screening sum rule for the underlying log-gas 
system, and is a fundamental requirement of S(9^\ 9 2 ) | Mar88 |. 

Substituting ( p. 18 ) in ( 2.16| ), interchanging the order of differentiation and integration ac- 
cording to ( 2.19) ), and using the Fourier expansion 



log|sin(^ 0) -9< 2 0) )/2\ 



00 



-27rlog2, 
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allows ( [2.16 ) to be evaluated as 



Var04)(°) = -5>a n c 



(2.20) 



n=l 



This has the well known feature of being independent of N — fluctuations are therefore strongly 
suppressed. Furthermore, it has been rigorously proved by Johansson [ J0I188 , Joh95] that the 
corresponding full distribution of the linear statistic A is given by the central limit-type theorem 

1 



Pi(u = A) 



(2.21) 



TV^oo (27T ( T 2 ) 1 /2 

where (A) is given by ( (2.15|) with the substitution ( 2.17 ), and a 1 = Yav(A) as given by ( 2.20 ). 

Let us now consider the situation for general fi, < 1. It is well known [Hua63|, and is 
simple to verify, that under the transformation 



e 3 — [x 



i0 



(0) 



(2.22) 



1 — fi*e J 

(note that the RHS has unit modulus and this mapping is one-to-one), the p.d.f. for general [i 
is transformed into the p.d.f. with fj, = according to 

N / 1 i 1 2 \/3(N-i)/2+l 



jjL*e ie 3 \ 2 



n 

l<j<k<N 

II le^-e^Vd^ 
l<j<k<N 



'ifd6x---de N 



deP. (2.23) 



This means that the correlation functions for general [i can be obtained from the correlation 
functions for fi = by applying the inverse of the transformation ( 2.22j ), 



(0) 



and noting that 



Hence, from ( 2.17 ), we have 

independent of j3, and so from ( 2.15| ) 

(A) = 



H + e i6 > 
1 + \i*e}®i 

(i-H 2 ) 



do. 



|1 - n*e ie \ 2 

n (1-H 2 ) 

2vr |1 - fi*e ie \ 2 



N 



2- 



2ir Jo |1 — fj,*e 



(i-H 2 ) 

* c i6\2 



p{6) dd. 



(2.24) 
(2.25) 
(2.26) 

(2.27) 



For the structure function, substituting ( 2.24| ) and ( |2.25 ) in ( 2.18] ) gives 



s(^ 0) ,4 0) )^i 0) 4 0) 



2 



log I e 



10 



(0) 



(o) 



1 



d 2 



(3tt 2 \d9td02 

1 / d 2 



Pit 2 V 00i 002 



log 



log 



J61 



J02 



1 + fi*e Wl \\l + fi*e id2 \ 

1^01^2. 



d9id6 2 



101 l6» 2 

e — e 



(2.28) 
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Thus S (61,62) for general fi is identical to 3(6^,6^) for /x = 0, and consequently 

Var(^) (0) = Var(A). (2.29) 

This fact illustrates a universality feature of the underlying log-gas: Var(^4) is invariant with 
respect to the particular one body potential modifying (LI), provided the corresponding one 
body density is a well behaved function. 



Finally, we note that (2.29) can be demonstrated via a numerical experiment. The ex- 



periment is performed by first generating random unitary matrices with uniform distribution 
(matrices from the CUE). This can be done by diagonalizing matrices from the GUE (random 
Hermitian matrices): the matrix of eigenvectors, when multiplied by a diagonal matrix with 
entries e lSj , (j = 1, . . . ,n) where 9j is a random angle between and 2-7T with uniform distri- 
bution, gives a matrix belonging to the CUE. Next we calculate the eigenvalues of each matrix 
(which will have distribution (LI) with j3 = 2), and transform them according to ( 2.22| ) with 



a specific value of [i. The resulting eigenvalues will have distribution as on the l.h.s. of djp . 
For each set k of eigenvalues {e iej }j=i t ... t N we then calculate := J2jLi a (@j) f° r some par- 
ticular choices of a. From the resulting list of values {A/%}, the empirical mean and standard 
derivation are calculated according to the usual formulas. In table 1 we present the result of 
performing this numerical experiment with N, the dimension of the unitary matrix, equal to 
15, and a(9) = cos jO (j = 1, . . . , 5). These empirical values are compared with the theoretical 
prediction for the variance in the limit N —* oo as given by (|2.2q ) (note that with a(9) = cos j6, 
a,j = ci-j = 1/2, a n = otherwise; thus ( |2.20|) gives Var(A) = j/2). 





Empirical 


Theoretical 


a(8) 


Vai(A) 


Var(A) 


cos 9 


0.509 


0.5 


cos 29 


0.972 


1 


cos 39 


1.6 


1.5 


cos 49 


1.8 


2 


cos 59 


2.6 


2.5 



Table 1. The second column contains the empirical variance of the quantity X^Li a (^j)> with 
a{9) as specified. This was calculated for 500 15 x 15 matrices with eigenvalue distribution given 
by the l.h.s. of ( [2.23 ) with (3 = 2 and /j, = .5. The final column contains the theoretical variance 
for the same quantity in the N — > oo limit. 



3 Laguerre ensemble with an TV-dependent exponent 
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3.1 Motivation 



Recently Brouwer et al. [BFB97] have considered the problem of the distribution of the eigen- 
values of the Wigner-Smith matrix Q = —ifrS~ 1 dS/dE. Here S refers to the scattering matrix 
coupled to a perfect lead which supports N channels of the same energy E. It was found that for 
each of the three possible symmetries of S, orthogonal {(5 = 1), unitary {(3 = 2) and symplectic 
{(3 = 4), the p.d.f. for the reciprocal of the eigenvalues of Q is given by ( |1.2| ) with a = N. 
This motivates a study of some of the properties of the distribution functions and fluctuation 
formulas associated with (|1.2|) for general a = YN, Y > 0. 



3.2 (3 = 1 and 2 

As remarked in the introduction, the p.d.f. (|1.2|) for (3 = 1, 2 and 4 is realized as the eigenvalue 
p.d.f. of random Wishart matrices A = X'X, where X has dimension M x N. For (3 = 1 
and 2, and with a proportional to N, the limiting form of the global density and the statistical 



properties of the largest and smallest eigenvalues have been extensively studied (see ||Ede88|| and 
references therein). In particular, with 

a = YN (3.1) 
it is known that the global eigenvalue density is given by 




where 



limp(4iVx)H -^^-^ W<*<W (3 . 2) 

N->oo 0, otherwise 



t 1 (Y)= 1 -(VTTY-lf, t 2 {Y) = \{Vl + Y + lf. (3.3) 



In qualitative terms, the result ( jjT^ ) says that the support of the density, which is (0, 1) 
when Y = 0, is repelled from the origin and elongated as Y increases. This is consistent with 



the log- gas interpretation of (|l.2j ) with a given by (3.1), as then there is an external charge of 
strength YN placed at the origin. This charge repels the N mobile charges of unit strength 
away from the origin. The fact that ([T^) is independent of (3 is also consistent with the log- 



gas interpretation. In fact macroscopic electrostatics says that the one-body potential in (1.5) 
results from a neutralizing background charge density Pb(y) according to 



x YN 
2 2~ 



logx + C = J^p b (y)log \x - y\ dy, x£l (3.4) 



where I is an interval in R + . The quantity (3 does not appear in this equation, so Pb{y) is 
independent of (3. But to leading order the particle density will equal the background density, 



as in general Coulomb systems strongly suppress charge fluctuations [Mar88|. The expected 
independence of (|3.2|) on (3 follows. 
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3.3 (3 even 



For (3 even, an exact /3-dimensional integral representation of the density in the finite system is 



available [ For94 , BF97a], which allows the global density limit to be taken explicitly. 

Now, in a system of N + 1 particles, the one-body density p(x) in the Laguerre ensemble is 
given by 

PN+i(x) := J** 1 e-W 2 xM 2 I N (a,P;x) (3.5) 

where 

N+l roo 

Z N+1 (a,P) := J] / dxie-^ 2 xf a/2 [J \x k - Xj f (3.6) 

1=1 l<j<fc<JV+l 
N roo 

I N (a,(3;x) := n/ dx^x-xtfe-^xf^ 2 \x k - Xj f (3.7) 

Z=l l<j<k<N 
The normalization (|3.6| ) is a well known limiting case of the Selberg integral, and has the exact 
evaluation (see e.g. [ Ask80[ |) 

_ , m ( 2 \ WHW-D/a y r(l + P(j + l)/2)r(l + + a)/2) ,„ 

ZAr(a5/?) = U n fxrr^ • (3 - 8) 



Notice that for (3 even the integral ( |3.7| ) is a polynomial in x. In this case ( |3.7D has been 
shown to be expressible in terms of a certain generalized Laguerre polynomial based on Jack 
polynomials |For94]. Furthermore, this generalized Laguerre polynomial has a different integral 
representation, which allows the ^-dimensional integral ( |3.7D to be expressed as a /3-dimensional 
integral. This reads |BF97a] 

' I N {a,/3;x) = —^—\f(a,f3;x)\ (3.9) 



Z N (a + 2 } (3) JVV ' Q(a,p) 
where 

f(a,(3;x) := 

I dti---[ dt l3 f[e xt nj N - 3+2/l3 (i-t j T+ N + 2 /f 3 - 1 n (tk-ttfW. (3.10) 

JCl Je P j=l l<j<k<0 

The contours of integration in ( |3.10[ ) must be simple loops which start at x = 1 and enclose 
the origin, and the quantity Q(a, (3) in (|3.9j ) is chosen so that at x = 0, the r.h.s. equals 1 (the 
l.h.s. has this property). Also, the modulus sign in (|3.9| ) has been included for convenience to 
eliminate terms of unit modulus which otherwise occur in Q(a,j3); this is valid for x G R since 
then the l.h.s. is positive. Choosing each Cj in ( |3.10| ) to be the unit circle gives 

Q(a, 0) = {2^ fMp{a + 2/(3 - 1, N, 2/(3) (3.11) 

where 

M N (a',b',c') := II / e me ^'^\l + e 2md r' +b ' R \e 2 ^ - e 2m9 ^ \ 2c ' 
" T(a' + b' + l + (j-l)c'ni + c'j) 

11 r(a' + 1 + (j - l)c')T(b' + l + l)d)T(l + d) 1 • ' 
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a=YN 



, which from (ED^) and 



(the integral is due to Morris [ Mor82f| ). 

Our interest is in the large- ./V asymptotic form of pn+i{ANx) 
( |3.9| ) requires the large-iV form of f(YN,(3;x). The necessary technique is a generalized saddle 
point analysis as introduced in |For94[| and presented in detail in [ BF97a |, The saddle points 
occur at the stationary points of 



NUxtj - log tj + (Y + 1) log(l - ij)) , 



(3.13) 



which is the ^-dependent term of the integrand of f(YN, (3\ x) when expressed as an exponential. 
A simple calculation gives that there are two stationary points, t+ and t_ say, given by 

x - Y/A ± ((x - Y/A) 2 - x) 1 ' 2 



Note that t+ = t*_ for 



2x 



(x - Y/A) 2 - x < 0. 



(3.14) 



(3.15) 



Assuming ( p. 15 ), following the strategy of [BF97a], the leading large- N asymptotic behaviour 
is obtained by deforming (3/2 of the contours through t + , and the remaining (3/2 of the contours 
through t- (this introduces a factor of (^2) *° accoun t f° r the number of ways of dividing the 
(3 contours into these two classes). The calculation now proceeds in a conventional way, with 
the exponent (3.13) being expanded about t± to second order, and the TV-independent terms in 
the integrand replaced by their value at t + or i_ as appropriate. After this step we have 



\f(YN,(3;x) 



( ^ 
\P/2 



( ^ 
V/V2 



\92(t+,Y)f 



dti • • • 



dt 



x n \t k 

i<i<fc</3/2 



|4//3 2 



-»g 1 (t + ,Y)(t*+-+t 2 fj/2 ) 



\92 



{t + ,Y)f\N gi {t + ,Y)\-^\v m {<. 



where 



9i(t+,Y) 

92(t+,Y) 
Vn(c) 



1 l + Y 

4"(i-M 2 

e 7V(4xt + -logt++(y+l)log(l-t + )) t -3+2//3^ 1 _ t+ ^/,8-l^ + 

/OO P 00 

dt x - dt N e-^l+-+t%) TJ \t k - tj \ 2c 



l<j<k<N 



N-l 



w n/2 n 



r(i + c(j + 1)) 
r(i + c ) 



(3.16) 

(3.17) 
(3.18) 

(3.19) 



The equality in ( 3.16| ) follow from a simple change of variables, while ( |3.19 ) is known as Mehta's 
integral, and can be evaluated as a limiting case of the Selberg integral [ Ask80 ] (for a direct 
evaluation see f|Eva94|l ). 

From the definition ( p. 14 ) we have that 



Ax' 



x+\ 



l + Y 

Ax 



x- 



^(x- (x-Y/Af 
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and use of these results in the formulas ( 3.1?1 ) and ( |3.18| ) defining g\ and 52 shows 

9 f +, lS-i = e 2 ^^™ y / 4 )(l+y)( 1+y )^/ 2 (4x)- y ^/ 2 (l+y) 1 /2l( x _( x _y / 4)2) 1/2 . ( 3 .20) 



\9i(U,Y)\ 
Furthermore, from 



Z N (a + 2,P) /2xiV/3 



Z N+1 (a,/3) \/3J 



r(f(a + JV + l) + l) 



r(i + f(i\r + i))r(fo + i)r(f(a + i) + i) 



'/9 



(3.21) 



while use of ( |3.12| ) in ( |3.11| ) and comparison with ( |3.19| ) shows that 
1 1 



Q(a,/3) (2ir)WVp(2/P) /J r(a + TV + |j) 



and use of ( [3.19D gives 



(V>/ 2 (2//3)) 2 _ ta /0 ^ /2 £(l+|) 



^(2//3) 



C9/2y 



r(i + /3) 



(3.22) 



(3.23) 



With ( |3.9|) substituted in (p75|) , the remaining task is to use Stirling's formula to compute 
the leading large-iV asymptotic behaviour of ( |3. 21 ) and ( |3.22| ) with a = YN. Performing this 
task, and substituting the resulting expression together with (|3~20| ) and (gj| ) in shows 
that for x such that ( p. 15 ) is true 

lim p N+1 (Wx) =-L( x -(x- y/4) 2 ) 1/2 . (3.24) 

N^oo a=YN Z7TX \ ' 

Outside the interval Q3.15D , i.e. outside x G [ x (l + ^ - vT+T), \(l + \ + vTTT)], this 
limit must vanish. This is seen from the fact that the density is positive, and must satisfy the 
normalization 

f°° 1 
/ p N+ i(ANx) dx ~ -, 
JO 4 

which is satisfied by the r.h.s. of ( p. 24 ). The result (|3.24j ) for the scaled density, established for 

even P, is identical to the result (p.2|) known for P = 1 and 2, as expected. 



4 The distribution functions in the neighbourhood of 
the smallest eigenvalue 

4.1 The n-point distributions 

For fixed N and p = 1, 2 or 4, the exact expressions for the general n-point distribution function 



in the Laguerre ensemble are known [ Bro65 , NW91 , NF95 |. With a = YN and N large, it is 



natural to move the origin to the (mean) location of the smallest eigenvalue, and to scale the 
eigenvalues so that the mean spacing near the spectrum edge is O(l). One anticipates that the 
limiting n-point distribution function will correspond to the n-point distribution function for the 
so-called soft edge, which is the edge of the spectrum for the Gaussian random matrix ensemble, 
with the eigenvalues appropriately scaled. 
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In quantitative terms, we expect that for appropriate v{N) independent of (3, 

lirn^ (^(iV)) n p {n) (iV(l - VT+Y) 2 - u(N)x h ...,N(l- VT+Y) 2 - u(N)x n ) 

= p$(x 1 ,...,x n ). (4.1) 



On the l.h.s., pr n \ refers to the n-point distribution function for the Laguerre ensemble with 
a = YN. Note from (|]|) and fl3.24|) that N(l - y/1 + Y) 2 is the location of the smallest 
eigenvalue (to leading order in N). On the r.h.s. 



where here pt n \ refers to the n-point distribution function for the Gaussian ensemble defined by 
the eigenvalue p.d.f 



i N 

1 rr _fl^2 



c 



n 

i=i 



Pxf/2 



I/ 9 



n \ x k-xjf, 



(4.3) 



and (2N) 1 ' 2 is the leading order location of the largest eigenvalue. In this section ( |4.1| ) will be 
explicitly verified for (3 = 2, with the quantity v(N) shown to be given by 



i/(JV) = iV 1/3 



2 1 / 3 {VT+Y-1) 



(y 2 -(2 + y)(VTTT-i) 2 ) 1 /3 



(4.4) 



Once v{N) has been determined, the validity of ( [4.1D for n = 1 can be established by 
matching the one-body density in the neighbourhood of the smallest eigenvalue implied by 
( 3.24] ) with the asymptotic behaviour [For92] 



P(i)V) 



(4.5) 



(this idea is motivated by a similar procedure used in |Nis96, KF97a[ ]). Now the result Q3.24 ) 
implies that for all (3 



P{1) (N(l-VTTY) 2 -v(N)x) 



N— >oo 
x — * — oo 



(i + y)V4 V (N) 



7T (1 - ^T+Y) 2 



N 



(4.6) 



valid for <C x -C N/v(N). Substituting Q4.5|) and fl4~q) in ( pTT| ) with n = 1, and using 
shows that (4.1) is satisfied provided 



1/4. 



2 1 I 2 {VTTy-i) 



{(2 + Y)(Y 2 - VT+Y-l) 2 ) 1 / 2 



(4.7) 



which is readily verified. 



In preparation for verifying (44.) for [3 = 2 and general n, we first recall some formulas 
particular to that coupling. For the Laguerre ensemble ( |l.2j ) we have [ Bro65| , 



p\n)(xi, ...,x n ) = det[P N {xj,x k )] jtk= 



l,...,m 



(4.8) 
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where with L^(x) denoting the Laguerre polynomial of degree n, 



r, I \ / \a/2 -(x+y)/2 L n( X ) L N -l(v) ~ L n(v) L N -l( X ) 

PN(x,y) := (xy) 1 e ^ c N , 

x-y 

^a/2 p -{x+ y ) 12 c N L% (x)y L a N \y) - L a N (y)xL a N '(x) 
N + a x — y 



(xy) c 



c N 



r(i + iv) 
r(o + n) 



(4.9) 



Furthermore, we know that | For93 | 

fffiixu ...,x n ) = det[K soft ( Xj , x k )} hk 

where, with Ai(x) denoting the Airy function, 

Ai(x)Ai'(y) - Ai(y)Ai'(x 



X,...,n 



K soit (x,y) :-- 



x-y 



Comparison of ( f4,10| ) and (|4.8j ) , and the fact that ( |4.1| ) is valid for n = 1 once 
shows that to verify (iA) it suffices to establish the asymptotic formula 



x NY/2 e - X /2 L Y N N {x] 



x=N(l-yT+Y) 2 -u{N)X N-^oo 



k N (Y)Ai(X) 



(4.10) 

(4.11) 
is established 

(4.12) 



with v{N) given by (4.4) and [i fixed. There is no need to specify k^(Y), as its value is 



uniquely determined by Q4.5| ) and Q4.6j). Indeed, substituting ( 4 . 1 2| ) in ( [4.9[ ), then substituting 
the resulting expression in fl4.8|) with n = 1, setting x = y and comparing with (iA) for x — > — oo 
shows that 

= - ' - igl. (4.13) 

V > (1-V1 + ^) 2 CN 

In table 2 we give the numerical value of the ratio of the l.h.s. to r.h.s. of ( [4.12] ) with k]y(Y) 
given by ( 4.1 3| ) for various values of X,Y,N. 





(0,1) 


(1,2) 


(-1,2) 


50 


1.0426 


1.0567 


0.9929 


60 


1.0402 


1.0556 


0.9938 


70 


1.0383 


1.0544 


0.9944 


80 


1.0367 


1.0533 


0.9949 


90 


1.0353 


1.0523 


0.9953 



Table 2. Numerical value of the ratio of the l.h.s. to r.h.s. of ( |4.12| ) for various values of N 
(column on left) and (X, Y) (top row). 



The asymptotic formula ( [4.12 ) can be derived by utilizing the fact that y 
satisfies the second order differential equation 



-x/2 x {a+l)/2 L a( x ) 



„ /27V + a + l 1 



y + 



2x 



+ 



a 



Ax 2 



y = o. 



(4.14) 
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Substituting a = YN, x = N(l - y/1 + Y) 2 - v{N)X, shows that for large N Q4.14Q reduces to 

1 /„, „ ^,\(i/(iV)) 3 



2(1 - ^/r+T) e 



Y 2 - (2 + y)(l - y/l + Y)'' 



N 



-xy = 0. 



(4.15) 



With v(N) given by ( [4.4| ), this equation reads y" — xy = and its unique solution which decays 
as X — > — oo is y = /cjv(F)Ai(X), thus establishing (|4.12[). 



4.2 Distribution of the smallest eigenvalue 

In principle, knowledge of the re-point distributions allows the calculation of other statistical 
quantities such as the distribution of the smallest eigenvalue. This together with the above 
results implies that, after appropriate change of origin and scale, the p.d.f. of the smallest 
eigenvalue in the Laguerre ensemble with a = YN equals, in the N — > oo limit, the p.d.f. for 



the largest (or equivalently smallest) eigenvalue in the Gaussian ensemble (4.3). For (5 = 2, this 
can be explicitly verified from the non-linear equations characterizing the respective p.d.f. 's due 
to Tracy and Widom [TW94a| , [TW94b|] . 



Let E{x) denote the probability that the interval (0, x) in the Laguerre ensemble (|l.2j) with 
(3 = 2 contains no eigenvalues, and let 

a(x) = x-^-logE(x). (4.16) 

ax 



Then it has been derived in [TW94b] that <r(x) satisfies the non- linear equation 



(xo") 2 = 4x(a'f + a 2 + (2a + AN-2x)aa' 

+(a 2 - 2ax - ANx + x 2 ){a') 2 - Aa{a') 2 . (4.17) 

Also, let E(x) denote the probability that there are no eigenvalues between (x, oo) in the (infinite 

dimensional) Gaussian ensemble ( |4.3| ) with coordinates as in ( |4~2D . Then it was derived in 
|TW94aH that the quantity 

R( x ) ■= ^L\ og E(x) (4.18) 

ax 

satisfies the non-linear equation 

{R") 2 + AR' ({R') 2 - xR' + R) = 0. (4.19) 



As the p.d.f. for the smallest (largest) eigenvalue is simply related to E(x) (E(x)) by dif- 
ferentiation, we see that the d.e. ( |4.17| ) characterizes the p.d.f. for the smallest eigenvalue in 
the Laguerre ensemble with (3 = 2, while ( 4.19| ) characterizes the limiting form of the p.d.f. of 



the eigenvalue at the spectrum edge of the Gaussian ensemble. (Of course boundary conditions 
must be specified; these follow from the small (large) x behaviour of the density.) The results 
of the previous subsection imply that 

lim e(n(VTTY - l) 2 - u(N)x) = E(x). (4.20) 

N— >oo V / a=YN 
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Indeed, a straightforward calculation using the explicit form p,4|) of v(N) shows that after 
changing variables x t—> N(\/l + Y — l) 2 — u(N)x in (4.17) and introducing the function 



5 {X) := U(N) a ( N{VYTY ~ 1)2 " U{N)X ) ' (421) 
the d.e. ( 4.17 ) reduces down to fl4.19| ) with R = a. This is precisely what is required by ( 4.20 ). 
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